A Study of the S=l/2 Alternating Chain using Multiprecision Methods 
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In this paper we present results for the ground state and low- lying excitations of the S = 1/2 alter- 
nating Heisenberg antiferromagnetic chain. Our more conventional techniques include perturbation 
theory about the dimer limit and numerical diagonalization of systems of up to 28 spins. A novel 
application of multiple precision numerical diagonalization allows us to determine analytical pertur- 
bation series to high order; the results found using this approach include ninth-order perturbation 
series for the ground state energy and one magnon gap, which were previously known only to third 
order. We also give the fifth-order dispersion relation and third-order exclusive neutron scattering 
structure factor for one-magnon modes and numerical and analytical binding energies of S = and 
S — 1 two-magnon bound states. 



PACS numbers: 75.10.Jm, 75.10.Hk, 75.30.Ds 
I. INTRODUCTION 

The alternating Heisenberg chain (AHC) is a simple 
quantum spin system that can be used to model the 
magnetic behavior of a wide range of materials; Table I 
gives some representative examples of alternating chains. 
This model is a straightforward generalization of the uni- 
form Heisenberg antiferromagnetic chain, which is the 
most widely studied quantum spin system. The uniform 
S = 1/2 chain has a gapless excitation spectrum with 
a known dispersion relation and a rather complicated 
ground state which is characterized by strong quantum 
fluctuations, making it highly unstable to perturbations. 

The alternating chain generalizes the uniform chain by 
alternating the spin-spin interaction between two values, 
Ji and J 2 . Since the alternating chain Hamiltonian is 
rotationally invariant with respect to spin, the total spin 
is a good quantum number, and the (antiferromagnetic) 
ground state is a spin singlet. The translational sym- 
metry of the uniform ID chain however is reduced by 
dimerization, and the resulting system has a gap to the 
first excited state, which has S = 1. This lowest excita- 
tion is part of a "one-magnon" triplet band. The alter- 
nating chain has a rather complicated spectrum of states 
at higher energies, including multimagnon continua and 
bound states. 

The alternating chain is of theoretical interest as a sim- 
ple ID isotropic quantum spin system with a gap, which 
presumably is qualitatively similar to other more com- 
plicated systems such as integer-spin chains and even-leg 
5 = 1/2 spin ladders. The approach to the uniform chain 
limit is also of interest as an example of critical behav- 
ior. Finally, the alternating chain is useful as a model for 
the application of new numerical techniques such as the 
multiprecision approach introduced here. 

The alternating chain Hamiltonian is realized in nature 
in many materials that have two important but struc- 
turally inequivalent superexchange paths that are spa- 



tially linked, so that a series of spin-spin interactions of 
strength J% — J% — J\ — J2... results. Examples of mate- 
rials of this type are (VO) 2 P 2 7 and Cu(N0 3 ) 2 • 2.5H 2 
and various aromatic free-radical compounds JjJ . 

Alternating chains may also arise as a result of the 
spin-Peierls effect. In the ID Heisenberg antiferromag- 
net a spatial dimerization of the ion positions along the 
chain gives alternating interaction strengths, which re- 
sults in a lowering of the magnetic ground state energy. 
There is a corresponding increase in the lattice energy 
(the phonon contribution) which dominates at large dis- 
tortions. In the combined magnetic-phonon system an 
equilibrium is reached at a spatial dimerization that min- 
imizes the ground state energy. This spontaneous dimer- 
ization is known as the spin-Peierls effect, and the re- 
sulting magnetic Hamiltonian is an alternating Heisen- 
berg chain. Examples of spin-Peierls alternating chains 
in nature are CuGe03 and a'-NaV 2 OV 

Much of the recent interest in alternating chains arises 
from the observation of a spin-Peierls effect in CuGeC>3 
||. Many experimental studies suggest an alternating 
chain interaction in CuGeC>3 (see for example |J and 
references cited therein), although interactions beyond 
nearest- neighbor are also thought to be important. The 
observation of a two-magnon continuum in CuGeC>3 with 
an onset close to 2 A || , where A is the magnon energy- 
gap at the zone-center, has motivated recent theoretical 
studies of the continuum and two-magnon bound states 
in the AHC. Added impetus has come from neutron scat- 
tering studies of (VO) 2 P 2 C-7 which show that this 
material is dominantly an alternating chain and provide 
evidence of a possible S = 1 two-magnon bound state. 

In this paper we present a detailed study of the AHC, 
including many new results for the ground state and low- 
lying excitations. We begin by introducing the model 
(Sec. II) and reviewing previous studies (Sec. III). Pertur- 
bation theory about the dimer limit, which we find to be 
particularly well suited to studying the AHC, is intro- 
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duced in Sec. IV. A and applied to the ground state en- 
ergy Eq , excitation gap E gap and one-magnon dispersion 
w{k) in Sec. IV. B. Sec. IV. C summarizes analytical pre- 
dictions for the critical behavior of the gap and ground 
state energy. Sec.V presents our numerical results for 
energies; Sec.V. A gives Lanczos results, and Sec.V.B in- 
troduces a new numerical method for abstracting ana- 
lytical perturbation series from high precision numerical 
results. In Sec.V.B we use this new approach to give se- 
ries (based on L — 20 diagonalization) to 0(a 9 ) for Eo, 
E gapi and the zone-boundary energy Ezb- Previously 
these series were known only to 0(a 3 ). We also used 
this multiple precision method to determine the series 
expansion of u>(k) to 0(a 5 ). These high-order formu- 
lae are accurate over a wide range of alternations and 
should prove useful to experimentalists. The critical be- 
havior is studied in Sec.V.C, and we present relationships 
between the derivatives of Eo and E gap at the critical 
point as well as comparing our results to the proposed 
scaling behavior of these quantities. Two-magnon bound 
states exist in the AHC, which is a convenient model 
for the study of this type of excitation. We discuss the 
binding mechanism and give second-order formulas for 
binding energies in Sec. VI. Since neutron scattering can 
give detailed information on the excitations of alternat- 
ing chains, we derive general expressions for the exclusive 
neutron scattering structure factor S{k) to a specific ex- 
citation (Sec. VII. A). We apply these results to the exci- 
tation of the one-magnon band in Sec. VII. B, and use the 
multiprecision method to calculate this S(k) to 0(a 3 ). 
A short discussion of the rather complicated neutron ex- 
citation of the two-magnon bound state band is given in 
Sec. VII. C. In Sec. VIII we present an illustrative applica- 
tion of our new formulae to a real material, the alternat- 
ing chain compound (VO)2P20y, for which single-crystal 
neutron scattering data is available. Finally we summ- 
marize our results and present our conclusions in Sec. IX. 



II. THE MODEL 



The AHC Hamiltonian is 

L/2 

H = Jl S%i-\ ■ §2i + J2 S21 ■ §2i+l 



(1) 



In this paper we impose periodic boundary conditions, 
with spins 1 and L + 1 identified. We usually assume 
that Ji > J 2 > 0, so we are in a regime of coupled 
antifcrromagnetic dimers. 

We can also write this in terms of J\ = J and the 
alternation a, where J2 = aJ; 



Nd is the number of independent dimers or unit cells, 
which are coupled by the interaction aJ. 

An equivalent form often used in the discussion of spin- 
Peierls transitions writes this as interactions of strength 
J{l+8) and J{\ — 6), which are related to our definitions 
by J = (1 + a) J/2 and S = (1 - a)/(l + a) . 

For a — 1 this system is an isotropic, uniform, S = 1/2 
Heisenberg chain which has gapless excitations, and for 
a = it reduces to uncoupled dimers with E gap = J. 
Since this is an isotropic Hamiltonian with antiferromag- 
netic couplings, for a > we expect an S = (singlet) 
ground state and an S = 1 (triplet) band of magnons as 
the first excitation. 

The geometry of our alternating chain is shown in 
Fig.l. Note that the unit cell has length b; this leads 
to a different set of momenta than the more familiar uni- 
form chain, which has a unit cell of a — b/2. Since the 
Hamiltonian is invariant under translations by multiples 
of b, the allowed momenta are 



k n — 



L/A b 



(3) 



N d =L/2 

H = J S2i~l ■ §2i + aJ £>2i ■ S; 



■ <->2i+l • 



(2) 



For L/2 even the index n takes the values n = 0, ±1, 
±2, . . ., ±(i/4 - 1), L/4; for L/2 odd the series stops 
with ±int(L/4 — 1). There are L/2 = Nd independent 
momenta because there are N4 invariant translations of 
H. Positive and negative k levels are degenerate as usual 
due to reflection symmetry. 



III. PREVIOUS STUDIES 

Early numerical studies of the zero temperature alter- 
nating chain by Duffy and Barr Q and Bonner and Blote 
considered the ground state energy and triplet gap 
on chains of up to 10 and 12 spins respectively. They 
concluded that this system probably had a gap for any 
nonzero alternation. Duffy and Barr also gave results for 
the ground-state nearest-neighbor correlation function, 
magnetization in an external field, and triplet dispersion 
relation uj(k). A principal concern of Bonner and Blote 
and subsequent numerical work was to test the critical 
behavior of the uniform chain limit; analytical studies 
had predicted that the gap E gap /J should open as <5 2 / 3 
times logarithmic corrections for small alternation [||, 
and that the bulk-limit ground-state energy per spin ex- 
pressed in terms of J and 5, e = E a /JL = 2e /(l + a), 
should approach 1/4 — ln(2) as S 4 ^ 3 times logarithmic 
corrections pLpOj. The dependence of eo on 6 is impor- 
tant in determining the existence of a spin-Peierls transi- 
tion in an antiferromagnetic chain coupled to the phonon 
field §. 

Numerical studies on larger systems were subsequently 
carried out by Soos et al. Jll]] (to L = 26 for eo and 
L = 21 for E gap ) and Spronken et al. £l|] (to L = 18). 
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Spronken et al. supported the anticipated critical behav- 
ior. Soos et al. however considered much smaller 5 and 
larger lattices and concluded that the expected asymp- 
totic form was incorrect. This issue is unresolved and 
merits future study on much larger systems. 

More recent studies of the alternating chain model were 
motivated by experimental work on CuGeC>3 In 
particular the question of possible two-magnon bound 
states has been of interest; an analytical paper by Uhrig 
and Schulz [|l3| anticipates an S = bound state for all 
6 and an S — 1 bound state "around k = 7r/2" (our k = 
7r/6, the zone boundary) "for not too small (5." . Bouzerar 
et al. jlj] similarly conclude that the S = 1 two-magnon 
bound state only exists for a range of k around the zone 
boundary. Fledderjohann and Gros Jl5| have searched for 
evidence of such bound states in a numerical study of the 
structure factor S(k,u}) on chains of up to L = 24, and 
conclude that an S = 1 two-magnon bound state does 
indeed lie below the two-magnon continuum for all S. 

Numerical studies of the thermodynamic properties 
of the alternating chain have received much less atten- 
tion. Duffy and Barr gave results for the internal energy, 
entropy, specific heat and magnetic susceptibility of an 
L = 10 chain for a range of alternations. Diederix et 
al. |l7j specialized to the parameter a = 0.27 appropriate 
for Cu(NC>3)2 • 2.5H2O and gave results for the magne- 
tization, susceptibility and entropy on systems of up to 
L = 12. Barnes and Riera (l8) gave results for the sus- 
ceptibility on chains of up to L = 16, and extrapolated 
to the bulk limit for values of a ss 0.6-0.8 considered 
appropriate for (VO)2P2C>7. 



IV. ANALYTICAL RESULTS 

A. Dimer Perturbation Theory 

Analytical results for the alternating chain can be de- 
rived using perturbation theory about the isolated dimer 
limit. For this purpose we partition the Hamiltonian into 
a dimer Hq and an interdimer interaction Hj, 



E 



J &2i- 



(4) 



The ground state of the full Hq is a direct product of 
S = dimer ground states, 



(0)\ 



N,, 

n 



= |o>, 



(6) 



with an energy of E = N d ■ E^ dimer) = -3JL/8. 

Similarly, the unperturbed one-magnon state with mo- 
mentum k and S z = +1 is given by 



i4 0) (fc)> (+) = -^E efc 



\{ + )m) 



(7) 



m— 1 



(We will suppress the redundant polarization superscript 
on I'i/'i) subsequently.) We take the location x m of dimer 
to to be the midpoint of the two spins. In this and similar 
state vectors, if the state of any dimer 77 is not specified 
explicitly it is in the ground state | o„). 

It is useful to derive the effect of Hi on dimer product 
basis states. For example, operating with Hi on the Hq 
ground state Eq.@ gives 



i/ / |0) = ^y3^|(0,0) m , m+1 ) 



(8) 



m — 1 



double excite 



where 



1(0,0) 



m ,m-|-l / 



1 

71 



!(+)*»(- Wi) - l(0)m(0) m+ l> + |(-) m (+)m+l) 



(9) 



is a state of two neighboring S = 1 dimer excitons at 
dimer sites 777, 777 + 1 coupled to give (S, S z ) = (0,0). 
Similarly the effect of Hi on a single S z = +1 exciton at 
dimer site m gives 

Hi \(+) m ) = y{ r l(+) m -i)-K+Wi> 

hop 



-|(+)m-l(0) m ) + |(0) TO _i(+) T 



excite 



Hi = 22 a J ■ §2i+l 



(5) 



-|(+)m(0)m+l) + |(0) m (+)m+l 



The single-dimer eigenstates of H are an S = ground 
state |o) = (|T|> - UT»/V2 with E^ mer) = -3 J/4 
and an S = 1 triplet of dimer excitations "excitons" 

{|(+)>. 1(0)), [(-)>}, with E[ dmer) = +J/4. We label 
these excitations by the dimer S z , for example |(0)) = 

(IU> + l!T»A/2. " 



N, 



+ v / 3 V'|(+) m (0,0) 



m' — l 



(10) 



double excite 



The prime on the sum indicates that all dimer sites 
represented in the state are distinct, so in this case 
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rn' ^ to, to — 1. Evidently Hj both translates the ex- 
citon (leading to momentum eigenstates) and couples it 
to two-exciton and three-exciton states of higher unper- 
turbed energy. The specific polarization state (|(+)(0)) — 
|(0)(+)))/v / 2 is forced because this is the unique \S — 
1, S z = 1) combination of two S = 1 dimers. We abbre- 
viate this state as |(l,l) mi „), specifying the Stotal and 
S z total and the excited dimers to and n (m <n), which 
gives the simplified form 



Hi \(+) m ) 



aJ 



|( + )m-l) + |(+)m+l> 



+V2 (j(l,l) 

rn — l,m 

> + 1(1,1) 

m,m+l/ I 



N d 



V3 ^'|(+) m (0,0) 

m' ,m' + l/ 



(ii) 



We can use this formalism to generate an expansion in 
a for the ground state and excitations and their matrix el- 
ements using standard quantum mechanical perturbation 
theory. These results are presented in the next section. 



B. Perturbative results for Eo, E gap and uj(k) 

The perturbative generalization of the ground state 
Eq.(i) to 0(a 2 ) is 



|V>o) = Va 



10) 



f F\ Nd 1 

^ ni—1 ' 



l(o,o) mi , 



l(+)m 1 (0) m2 (-) m3 ) + |(0) mi (-) m2 ( + )„ l3 



+l(-)mx(+)m 2 (0)m 3 ) - I (+)m t (-)m, (O)m,) 



-K-k(oU(+) m 3)-l(oW(+U(-U)J , (is) 

is the unique S — combination of three spin-one objects 
at adjacent sites. 

Since the 0(a p ) state determines the 0(a 2p+1 ) energy, 
we can in principle use Eq. (|l2| ) to derive the ground state 
energy to 0(a 5 ). This proves to be a rather intricate 
calculation. We have carried out this derivation of eo = 
Eq/LJ analytically to 0(a 4 ), with the result 



e (a) 



3 3 9 3 q 13 a 

a a a 

2 3 2 6 2 8 2 12 



(14) 



This series was previously evaluated to 0(a 3 ) by Brooks 
Harris g§. 

A similar O(a) generalization of the unperturbed 
S = 1 one-magnon excitation Eq.(^) gives 



\Mk))=Vi 



r i ^ 

v v m— 1 



//•:/) 



rn— 1 



1) e l 



Akx n 



\( + U 



(1,1 



j J- Jm, 771+1/ 



/o N d /o N d 

x El(o.oU + i)^El(o.oW^) 



32 



m— 1 



i 



771 — 1 



16 V 3 



1(0,0)™,™+!, 



rn+2/ 



rn— 1 



N,, 



128 S '!<».») 



m,m+l 



(0,0) 



m' ,m' + l/ 



(12) 



where 770 = 1 — (3 / 128) a 2 Nd is the 0(a 2 ) normalization. 
Note that three- and four-exciton states appear at 0(a 2 ). 
The four-exciton states encountered here are two S = 0, 
(0, 0) mjm +i two-exciton pairs, again with the restriction 
on the sum J2 that no excitons overlap. The thrcc- 
exciton state with (S, S z ) = (0, 0), 



V3 



\ ^ I Akx 



8 *-f 

m.rn — 1 



l(+)m(0,0) r 



(15) 



The normalization r)i = 1 to this order in a. Taking the 
expected value of H with this state gives a one-magnon 
dispersion relation of 



w(fc) 



1 1 2 3 
1— — a +— a 



1G 



64 



~ ( ili " 2 + ^ " 3 ) cos ( 2fc6 ) - 64 " 3 cos ( 3fc6 ) 
The (9(a 3 ) gap is therefore 



E„ 



J 



I- -a- 
2 



32 



■ a 



(16) 



(17) 



These one-magnon energies were found previously to this 
order by Brooks Harris |l9), and serve as a check of our 
0(a) one-magnon state Eq.(|T5|). 
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C. Critical behavior 



which implies for any energy eigenvalue 



As we approach the uniform chain the ground state 
energy and the gap are both expected to approach their 
limiting values as powers of 6 times logarithmic correc- 
tions |{L(J. The behavior near the critical point is usu- 
ally discussed in terms of the variable 



with 



5=(l-a)/(l + a) , 



J = (1 + a) J/2 



(18) 



fixed, so the alternating couplings are J(l + 6) and 
J(l — 5). These variables are more appropriate for a 
spin-Peierls system because a displacement of an inter- 
mediate ion by 0(S) should increase and decrease al- 
ternate couplings by approximately the same amount. 
The ground state energy per spin relative to fixed J is 
e (<5) = 2eo/(l + a). 

The critical behavior of the ground state energy and 
singlet-triplet gap has been discussed by Cross and Fisher 
H and by Black and Emery . The approach used was 
to consider the properies of the Heisenberg chain within 
a Luttinger-Tomonaga approximation, which involves a 
Jordan- Wigner transformation to a fermion representa- 
tion of the spin operators, and then replacing the cosi- 
nusoidal fermion dispersion by a linear dispersion at the 
Fermi wavevector. This linear approximation is required 
to simplify the commutation relations between the den- 
sity operators, allowing the interacting fermion problem 
to be solved. Renormalization techniques are then used 
to calculate the asymptotic behavior of various physi- 
cal quantities within this approximation. The approach 
makes uncontrolled approximations by neglecting states 
far from the Fermi surface and ignoring energy renor- 
malization effects, so it is unclear whether this gives the 
correct critical behavior. The predicted asymptotic 5 de- 
pendence is 



lim en (S) — en(<5 

5^0 



0) oc 



llnJI 



and 



lim^ 

5^0 J 



OC 



ln<5|V2 



(20) 



(21) 



We will compare these predictions with our numerical 
results in the next section. 

One may derive some relations between energies and 
their derivatives near the critical point from a simple 
identity satisfied by the alternating chain Hamiltonian 
Eq.(B). Note the proportionality relation 



H(J,aJ) = a ■ H(a~ 1 J, J) 



(22) 



E n {a) 
J 



E n (a 



J 



(23) 



Assuming that there are no singularities on the real 
axis except at a = 1, we may differentiate this relation 
with respect to a elsewhere. As it is expected that the 
singularity in eo(a) = Eg/ JL at the critical point a = 1 
is higher order than linear, deo(a)/da should be well 
defined everywhere; differentiating Eq.(23) with n = 
therefore leads to 



(19) de (a) 



da 



1 , 
= 2 e °( a 



1_M2) 

> O 1 



= -0.22157.. 



(24) 



This is consistent with the expectation that the scaled 
en(<5) has zero slope in 5 as we approach the critical point. 
To see this, note that Sq(S) = 2en(a)/(l + a), so 



de (S) 



e (a) 



cfen(a) 



d5 " ' (1+5) da 

and as we approach the critical point 



,. de (S) 
hm — — - 

<5->o do 



lim en (a) 



de Q {a) 
da 







(25) 



(26) 



Successive derivatives of Eq. (|23|) can be used to infer rela- 
tions between higher derivatives of en (a) (or other energy 
eigenvalues) as one approaches a = 1. 



V. NUMERICAL METHODS 
A. Lanczos results 

Direct numerical diagonalization of moderately large 
systems is possible for the 5 = 1/2 alternating chain. 
Here we used a Lanczos method |2^] to obtain ground 
state and one-magnon energies on L = An lattices up to 
L = 28. Motivated by previous numerical studies, we 
extrapolate these energies to bulk limits using a simple 
exponential-and-power estimate for the finite size depen- 
dence, 



f(a,L) = f(a) + ci 



exp(~L/c 2 ) 
LP 



(27) 



where p = 1 for energy gaps (/ = E n — Eg) and p = 2 
for the ground state energy per spin (/ = Eq/L). Wc 
determined the finite-lattice energies to about 14 place 
accuracy, and fitted the L — 4(n — 2),4(n — 1) and 4n 
results to these asymptotic forms. 

The resulting bulk-limit ground state energy is shown 
in Fig. 2 and presented in Table II, to nine figure ac- 
curacy for the smaller a values. For larger a we in- 
clude the change in e (a) between L — 16,20,24 and 
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L = 20, 24, 28 extrapolations in parenthesis after the tab- 
ulated L = 20, 24, 28 result, as an error estimate. These 
numerical energies provide an accurate check of the per- 
turbative formulas Eq.(|l4|) and Eq.(p8|). 

We also used Lanczos diagonalization on lattices up 
to L = 28 to determine the singlet-triplet gap and zone 
boundary energy. These results are given in Table II and 
shown in Fig. 3, again with a systematic error estimate 
that is the discrepancy between L < 24 and L < 28 ex- 
trapolations. (The + sign indicates that the bulk-limit 
gap estimate increased with increasing L.) Our Lanc- 
zos results are again consistent with the perturbative ex- 
pansions Eq.([l6|) and Eq. (|l7|) . The higher-order multi- 
ple precision series Eq.(|29|) and Eq. (|30|) are as expected 
found to be in agreement to much higher accuracy. 



B. High-order series from multiple precision 

One may use numerical diagonalization combined with 
multiple precision programming to determine analytical 
perturbation series to high order. For this novel ap- 
plication of numerical methods to spin systems we em- 
ployed the multiprecision package MPFUN developed by 
D.Bailey ]2~H ], applied to our Fortran code for low-lying 
eigenvectors of the alternating chain using the "modified 
Lanczos" method We typically generated energies 

to 300 significant figures with a — 10~ 30 (and to 400 fig- 
ures for L = 20), which allowed the perturbation expan- 
sion coefficients to be read directly from the numerical 
energies as rational fractions. This was possible in part 
because the energy denominators involved simple powers 
of small integers that could be anticipated. The limiting 
order in this approach is determined by the size of the 
system one can diagonalize, since the periodic boundary 
conditions introduce 0(a L / 2 ) finite-lattice corrections to 
energies. This gave a limit of 0(a 9 ) for the order of the 
bulk limit expansion that could be determined from the 
largest system we diagonalized with multiple precision, 
L = 20. 

The multiprecision package was implemented on a Pen- 
tium PC, a DEC Alpha and a Sun 450. Execution times 
for Eq at this level of precision were approximately 6 CPU 
hours for L = 16 on the DEC Alpha and 100 CPU hours 
for L = 20 on the Sun 450. 

The 0(a 9 ) series for the ground state energy per spin 
determined in this manner is 



e (a) 



_3_ 



2 6 



a 



28' 



13 



2 14 -3 



463 

117 . I 



7-61-191 7 11-139-271 

-a — 



2 22 • 3 3 



107-22005559 
2 30 -3 5 -5 2 ' 



2 21 • 3 4 • 5 



(28) 



and the 0(a 9 ) series for the gap to the k = one-magnon 
state is 



j 2 « 2 3 " + 2 s " 2? . 3 



a 



761 



(ll) 2 -157 
2 16 • 3 3 



21739 
2 18 • 3 3 



107 ■ 283 ■ 707 9 
2 24 • 3 4 • 5 



1307-9151183 
2 28 • 3 6 • 5 2 



(29) 



The zone-boundary (k = n/b) energy of the one-magnon 
state relative to Eq to this order is 



E ZB n , 1 , 1 2 1 

= 1 H — a -\ — -a = 

J 2 2 3 2 5 • 3 



S3 



71 • 149 
2 16 • 3 3 ' 



6373 
2 14 • 3 4 



41 ■ 256687901 
2 28 • 3 6 • 5 2 ' 



19 • 128461 
2 24 • 3 2 • 5 



a 



(30) 



We can also use multiprecision methods to determine 
the one-magnon dispersion relation w(fc), parametrized 
by 



Lo(k) 



ae(a) cos(£kb) . 



(31) 



On a finite lattice, momenta are only defined at the Nd — 
L/2 independent values of Eq.(||). (See also Fig. lb.) The 
relation uj(k) = Lu(—k) further reduces this to a total of 
int(Nd/2) + 1 independent lattice energies. These can be 
expanded as power series in a, and we again encounter 
finite lattice artifacts in these expansions beginning at 
0(a N *). 

The use of the lattice {ui(k n )} to determine the Fourier 
coefficients in Eq. (pl|) is nontrivial because there are in- 
finitely many coefficients but only int(Nd/2) + 1 lat- 
tice energies. We can proceed by determining the 
int(Nd/2) + 1 finite lattice Fourier coefficients {at(a, L)}, 
I = 0, 1, ...,int(N d /2), defined by 



oj(k n 



int(N d /2) 

^ ae(a, L) cos(£k n b) 



(32) 



where k n b = 2%n/Nd- 

We can invert this using a general result for the sum 
of a product of cosines over the discrete lattice momenta. 
We shall now assume Nd is even, so there are singly de- 
generate k = and k = n/b points in addition to the 
doubly degenerate values k = ±2tt /Ngb, ±4ir/Ndb, .... 
We can translate the k = and negative k values by 
2ir/b, so a sum over lattice k values becomes a sum over 
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n = 1, 2, N d , with k n = 2irn/Ndb. The summed prod- 
uct of cosines is 

N d 



— ^ cos(2Tr£n/N d ) co8(2w£n'/N d ) 



2 (^mad(n—n',N d ),0 + $mod(n+ri ,N d ),oj ■ 



(33) 



On multiplying Eq.(|32|) by cos(2ttii£' /N d ) and summing 
over n, we therefore find the lattice Fourier coefficients 



For the other Fourier modes the order of perturba- 
tion theory to which the bulk-limit coefficients are de- 
termined is reduced by the contributions of higher-order 
coefficients. For £' = 1 for example we find 



ai 



2 to(k n ) 

a n=\ 



J 



cos{2im/N d ) 



a i (^mod(t-l,N d ),0 + dmod(i+l,N d ),o) 



ae(a,L) = 



«1 + a N d + l + 0-2N d + l + ■■■ + a N d -l + 0-2N d -l + 



(37) 



N d 



N d (1 + 5lo + k,N d /2) ~ 



(34) 



If we assume more generally that the {cu(k n )} are sam- 
plings of the continuous function given by Eq.(3_l), and 
invert the Fourier expansion using Eq.(^) we find 



A',; 



N d ^ J 

n— 1 



oo 



(35) 



i=0 



As anticipated this constitutes int(N d /2) + 1 constraints 
(the number of independent choices for £') on the infinite 
set of Fourier coefficients. 

As an illustration, for £' = Eqs.(pJ) and (|35|) imply 



a 



AT., 

n=l 



The contributions of the two delta functions that com- 
bined in the previous £' = example are now distinct, 
so there is an 0(a Nd ~ 1 ) difference due to aN d -i between 
the known lattice series for Si and the bulk-limit a\ se- 
ries. The worst case is to' = N d /2 — 1; the lattice and 
bulk-limit coefficients are then related by 

&N d /2-l = 

a N d /2-l + a ZN d /2-l + ••• + a N d /2+l + a 3AV2+l + ••• i 

(38) 

so the series for the bulk-limit a Nd / 2 -i, which begins 
at 0(a Nd / 2 ~ 1 ), cannot be determined beyond 0(a Nd ^ 2 ) 
due to the presence of a Nd / 2 +\- Thus we conclude 
that the perturbative expansion of the bulk-limit dis- 
persion relation Eq.(|3l|) can only be completely deter- 
mined to 0(a Nd / 2 ) from finite lattice data. The at- 
tainable order depends on the mode, and is in general 
O^-i-modj^/a)) for a e with N d /2 >£>0. 

It follows that the complete set of bulk-limit Fourier 
coefficients {at} is uniquely determined just to 0(a 5 ) 
by the diagonalization of the L — 20 lattice. These co- 
efficients, which confirm and continue the one-magnon 
dispersion relation of Brooks Harris, Eq.(|lrj), are 



2 5>(' 



'mod(t,N d )fi 



id(e,N, 



ao = 1 



«1 



16 " 



64 " 



^_a 4 

1024 



"256 a 



256 



2048 



ao + a-N d + a-2N d + 



(36) 



a 2 



512 " 



283 5 
18432 " 



Since we know that the perturbative series for both 
the lattice {de} and bulk-limit {at} coefficients begin at 
O(cZ), the expansions of ao and do must be identical 
until we encounter contributions from higher-order coef- 
ficients; for ao these begin at 0(a Nd ) due to a,N d . This is 
also the order at which finite-lattice artifacts appear in 
the energies, so there has been no loss of order reached 
in the ao expansion beyond the usual finite-lattice limi- 
tation. An equivalent conclusion follows for the Fourier 
coefficient a^ d / 2 . 



«3 



«4 



64 " 



-J- a 4 ^-a 5 

48 1024 



1024 



9216 



4096 



(39) 

The fifth-order one-magnon dispersion relation given 
by Eqs.(|3l],|39]) is shown for a range of alternations in 
Fig.4. 
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C. Numerical results for critical behavior 

The expected behavior of E gap /J and eo(5) = 
Eq/ JL = 2ep(a)/(l + a) as we approach the critical 
point S = is MM 



lim en (S) 

6^0 



eo(5 = 0) oc 



£ 4 / 3 
UnJI 



and 



lim^Zi 

s^o J 



§V3 



x 



(40) 



(41) 



We will now compare these theoretical asymptotic forms 
with our numerical results for Egap/J and eo(S) at small 



interest in two-magnon bound states in dimerized quan- 
tum spin systems. The existence of such bound states 
in the alternating chain model appears very likely, since 
the interdimer interaction Hj, Eq.(|^), gives an attractive 
diagonal potential energy between two adjacent excited 
dimers if their total spin is 5 = or 5 = 1. To see this, 
note the effect of Hj on a state of adjacent (+)(0) dimcr 
cxcitons; 



Hj\ o (+)(0)o) = 



aJ 



o(0)(+)o) 



transpose 



-|(+)°(0)o)-|o(+)o(0)> 



hop 



Taking the ground state energy per spin first, in Fig. 5a 
we show a logarithmic plot of en versus 5; evidently the 
smaller S values do support approximate power-law be- 
havior. A fit of the three smallest 6 points to the form 
c\8 p gives ci = 0.3632 and p = 1.412, apparently con- 
sistent with 5 4 / 3 times logarithmic corrections. The an- 
ticipated theoretical form ci<5 4//3 /| ln<5 | however does not 
give an especially good fit, and is actually worse than a 
pure 4/3 power law. A 4/3 power fit to the three smallest 
6 points gives c\ = 0.3134, which is reasonably accurate 
for larger 5 as well; this fit is shown in Fig. 5a. 

The 6 dependence we observe for the singlet-triplet 
gap scaled by J is shown in Fig. 5b. The approximate 
linearity of the log-log plot suggests power-law behavior 
over the full range of S considered. A fit of all points 
to ci8 p gives ci = 1.999 and p = 0.7497, suggesting 
a power law of <5 3 / 4 rather than the theoretical <5 2 / 3 of 
Eq. (^H) . Fitting ci(5 3 / 4 to the three smallest 6 points 
gives ci = 2.003, which is shown in Fig. 5b. Evidently 
this is a remarkably good fit. The surprising accuracy 
of this form for large alternation can be understood by 
noting that the similar function Egap/J = 2 5 3 / 4 , cor- 
responding to Egap/J = (1 — a) 3 / 4 (l + a) 1 / 4 , gives the 
correct E gap /J power series Eq.(|i"7j) to 0(a 2 ). 

The form c\S 2 l' i j\ In S I 1 / 2 found by Black and Emery 
gives a much less accurate description of our data. Thus 
we appear to support a different asymptotic power law 
than expected for the singlet-triplet gap, although we 
cannot continue to very small 6 to see if this discrepancy 
persists. The behavior of the ground state energy and 
gap much closer to the critical point would be an inter- 
esting topic for a detailed investigation on large lattices, 
using for example the DMRG approach pa]. 



VI. TWO-MAGNON BOUND STATES 

Recent theoretical work [pL 3| — [L 6|| motivated by neutron 
scattering studies of CuGeC>3 has led to considerable 



o (+) o o) + | oo(+)o) 



dccxcitc 



-|(+)(0)(0)o)-|(0)(+)(0)o) 



+1° (+)(+)(-)> -1° (+)(-)(+)> 



N d 

\/3 | (0, 0) mW +i •••(+) (0).. 

' * ^ 

m' = l 



(42) 



double excite 



Again in the index ml takes on all values that do not 
superimpose a dimer excitation on the already excited 
(+) or (— ) sites. 

The analogous expressions for other polarization states 
may be combined to give the effect of Hi on nearest 
neighbor dimer excitons with definite total spin. There 
is a diagonal interaction due to the "transpose" matrix 
element and others that retain two neighboring excitons. 
This static inter-exciton potential is 



\H!\(S,S Z ) 



m,m+l/ 



+Ja/A S = 2- 
-Ja/4 5=1; 
-Ja/2 5 = 0. 



(43) 

This suggests that at large alternation (small a) we 
should find two-magnon bound states with 5 = and 
5 = 1, with binding energies of approximately Ja/2 and 
Ja/4 respectively. (This only applies to k = 7r/fo; at 
other k values there is a complication that modifies this 
result, as noted below.) 

We can again use numerical results with small coupling 
a to establish the higher-order perturbation series for 
properties of these bound states. Unfortunately this is a 
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much more difficult numerical problem than the study of 
the ground state and one-magnon levels, so here we give 
only a few preliminary results. 

The binding energies Eb of the 5 = and 5 = 1 bound 
states are defined by 



EbQc) — min 

k±+k 2 =k 



(w(*i)+a;(fe)) - (E(k)-E 



(44) 



(The u> sum gives the onset of the two-magnon continuum 
at k.) Assuming that the continuum onset at k = n/b 
is given by k\ — and &2 = n/b, we find that the two- 
magnon binding energies to 0(a 2 ) on L > 12 lattices 
are 



E B (k = n/b) 



32" 



5 = 1, 



^a 2 5 = 0. 



(45) 



Our numerical result for the bulk-limit binding energies 
at k = n/b (Table II and Fig. 6) are clearly consistent with 
these perturbative results at small a. It is interesting 
that the binding is much weaker near k = 0, for example 
the expansion of the 5 = binding energy at k = 
appears to begin at 0(a 2 ). 

The determination of these perturbative binding ener- 
gies analytically at general k is a complicated problem 
because the two-exciton sector is a manifold of degener- 
ate states under H . Determining the appropriate basis 
states within this degenerate manifold requires diagonal- 
ization of the "hopping" part of Hj. This is only straight- 
forward at k — n/b, where this hopping amplitude van- 
ishes; as an illustration, for adjacent (+)(0) excitons we 
find 

( (+) o (0)o ; k\Hj\ o (+)(0)o ; k) cx cos(fc&/ 2 ) . (46) 

Thus k — n jb bound states do not mix with the two- 
exciton continuum to leading order, and remain relatively 
localized. 

Except near k = n/b we expect the coupling be- 
tween nearest-neighbor and separated excitons to be very 
important, because the energy denominator is 0(aJ), 
which is the same order as the hopping matrix element. 
This implies that there are corrections to the bound-state 
wavefunction of order 0((Hj)/A£) = 0{aJ)/0(aJ) = 
0(a) and corrections to the static binding energy 
Eq.@ of order 0(\ (Hi) \ 2 / AE) = 0(a) except at k = 
n/b. 

At the extreme point k = we see no evidence for 
a 5 = 1 bound state in our numerical extrapolation of 
Lanczos data. Although we do see evidence of an 5 = 
bound state at k = 0, it is quite weakly bound and (in 
consequence) has a very extended spatial wavefunction. 
The difference between k = and k = n/b bound state 
wavefunctions was quite evident in the finite-size effects 
seen in our numerical extrapolation. 



In a preliminary numerical study of a truncated sys- 
tem of zero-, one- and two-exciton states on an L = 200 
lattice we find that the attractive potentials in Eq.(|4^) 
are strong enough to form an 5 = bound state for all 
k, but apparently the 5=1 bound state exists only for 
a range of k around n/b. Similar conclusions have been 
reported by Uhrig and Schulz Jl3[ and Bouzcrar et al. 

& 

These results are especially relevant to the alternating- 
chain material (VO)2P2C>7, "VOPO". An excitation has 
been observed in VOPO just below the two-magnon con- 
tinuum, which has been cited as a possible two-magnon 
bound state [0,f|. Since this peak is seen clearly for a 
range of k including k « 2n/b (equivalent to the k = 
point where we find no 5 = 1 bound state), it is in- 
consistent with the expectations of the alternating chain 
model for such a bound state. If mode observed in VOPO 
is indeed a two-magnon bound state, its persistence to 
k In /b is presumably due to additional interactions. 



VII. NEUTRON SCATTERING STRUCTURE 
FACTOR 

A. General results 

Identification of the magnetic excitations predicted by 
the alternating chain model will be facilitated by esti- 
mates of their couplings to external probes, such as pho- 
tons (especially for 5 = states, through Raman scat- 
tering) and neutrons (for 5 = 1 states). Here we present 
perturbative and numerical results for the neutron scat- 
tering structure factor to the first two 5=1 excitations 
in the alternating chain, which are the one-magnon mode 
and the 5 = 1 two-magnon bound state. 

The neutron scattering cross section is proportional to 
the structure factor, which in the Heisenberg picture is 



,(k,w) 



1 

27 



dt e lujt e lk ' {i 



^2 (V'olS'mC^i)*)^'^) )!^) 



(47) 



sites i,j 



We can insert a complete set of eigenstates {\iI>n)} of 
the full H between the spins in this matrix element and 
write the spin structure factor as a sum over exclusive 
structure factors, 



N 



s-(^)=E s -"(^) 

where 



(48) 



sites i,j 



J2 ~ dt ^ e iHSi ~ Sj) 
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(lPo\ S L (% . I V'iv) (l/>N | 5m' (f i , 0) | ^o) 



(49) 



Each exclusive structure factor S*!^ (k, to), gives the 
intensity of scattering from |^>o) to a specific excited state 
\iPn), 

The Heisenberg picture operator S^Xj, t) — exp(iiii) 
S^XjjO) exp(-iHt) gives trivial exponentials in t, so 
the time integral leads to a single energy-conserving delta 
function. In isotropic antiferromagnets the ground state 
typically has 5 = 0, so the accessible excited states 
(IV'jv)} ah have 5 = 1. We can then evaluate this ex- 
clusive structure factor to a specific polarization state, 
here S z = +1, without loss of generality. For S z = +1 
the only nonzero spherical components of the exclusive 
structure factor is then 



5(E N -Eq-u) 



E 



(ip N \S + (xi)\ipo) e 



ik-. 



(50) 



If the excited state \i/jn) is an eigenstate of momentum, 
which we can assume for the alternating chain without 
loss of generality, the matrix elements of the spin opera- 
tor at translationally equivalent sites are equal modulo a 
plane wave, 



(ip N (p)\S m (xi)\ip ) 



-ip-(xi-So) 



(lpN(p)\S m (x )\lp ) 



(51) 



where xo is some reference site. The sum over all spin 
sites i can then be reduced to a sum over all sites in the 
unit cell i* times a momentum conserving delta function, 



Nunit cells $(En 



(^n(p)\ 5+(f;*) \ip ) e 



(52) 



If we are only interested in the k dependence and rel- 
ative intensities of neutron scattering from the various 
5 = 1 excitations, the overall normalization is irrelevant, 
and we can simply evaluate dimensionless reduced inten- 
sities. It is convenient to normalize this reduced exclusive 
structure factor as 



B. One-magnon exclusive S(k) 

We will now consider the reduced exclusive neutron 
scattering structure factor S(k), Eq. (|53|) , for the excita- 
tion of one-magnon states in the alternating chain, using 
dimer perturbation theory and the multiprecision tech- 
nique. The exclusive structure factor involves the matrix 
element 



(^N(k)\S+(x„) |Vo> e** 2 < 



-ikd/2 



5*e 



(54) 



The superscripts L and R refer to the left and right spins 
in the first dimer. Note that this can be written as a sum 
of dimer-spin conserving and changing terms, 

gL e -ikd/2 _|_ gR e +ikd/2 _ 



cos(fcd/2) (s$ + S$\ 



Iimer-spm conserving 



sin(fcd/2) 7 ( 5^ - 5* 



(55) 



dimer— spin changing 



For the one-magnon excitation it suffices to determine 
the matrix element of the raising operator on a single 
spin, because Si and 5+ have opposite one-magnon ma- 
trix elements for any k. (This is not true for a general 
5 = 1 excitation.) Evaluating the matrix element of Si 
to 0(a) using the perturbative \ipo) and \ipi(k)), Eq.(p^) 
and Eq.dl5|), we find 



(^l(fe)IS^I^) =-^= U + jC08(kb) 

so the matrix element Eq.(HJ) is 

]T (iP N (k)\S + (x„) |^o) = 



(56) 



S(k) = N, 



unit cells 



(ip N (k)\S + (x„) |Vo) e***" 



(53) 



(Here and for the remainder of the paper we suppress the 
superscript ipo — > ippj on the exclusive structure factor.) 
This expression can be evaluated analytically (using per- 
turbation theory or other approximate wavefunctions) or 
numerically using wavefunctions on finite lattices. 



e ifcd/2 _ g-ifcd/2 



1 H cos(kb) 

4 



(57) 



The one-magnon exclusive neutron scattering structure 
factor S(k) is proportional to the modulus squared of 
this spin matrix element. To 0(a) we find 



S(k) = ( 1 - cos(kd) ) ( 1 + - cos{kb) 



(58) 
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The small-/c suppression 1— cos(fcd) oc sin(fc<i/2) 2 is famil- 
iar from isolated dimer problems and gives a basic inten- 
sity "envelope" that measures the dimer size d. (See for 
example Ref. pifl .) The separation between dimer centers 
b enters as a more rapid modulation 1 + (a/2) cos(kb) of 
the intrinsic dimer form 1 — cos(fcd). This 0(a) mod- 
ulation arises from the excitation of the "two-exciton" 
component of the ground state \ipa) to the unperturbed 
"one-exciton" component of \ipi(k)), and apparently has 
been observed in recent neutron scattering experiments 
on single crystals of Sri 4 Cu2404i [ p5[ . 

The S(k) series Eq.(|58|) may be continued using the 
multiprecision approach. We introduce a cosine expan- 
sion for the Si matrix element, 



excited by a rather different mechanism than the one- 
magnon state discussed in Sec. VIII. B above. To leading 
order the bound state wavefunction is 



1(1,1, 



m,m+l/ 



(62) 



which is invisible to neutron excitation of the "bare" 
ground state, since two spin flips are required to connect 
these states. 

At 0(a) a coupling to neutrons appears through the 
nonleading parts of the ground state and bound state. 
There are two such contributions, which are apparent on 
inspecting the 0(a) states, 



(Mk)\S^ |^o> = f>(«) cos(£kb) , 



(59) 



\^o) = \0)-\a £ |(0,0) m , m+1 ) 

m— 1 



(63) 



and the coefficients we find to 0(a 3 ) are 



and 



s = 1 


11 

64 


a 2 


S- a 

128 


«1 = 


^4 " 16 


a 2 


+-^-a 

~ 1536 u 


S2 = 


^64 


a 2 


^384 " 


S3 = 






+^a 

^512" 



(60) 



These have the same order finite size effects as the dis- 
persion Fourier coefficients {at}, so diagonalization of an 
L site chain gives the complete set of exclusive structure 
factor Fourier coefficients to 0(a L / 4 ). 

These coefficients determine the 0(a 3 ) generalization 
of Eq. (|58|) , which is 



S(k) = 1 - cos(fcd) 



'l_JL a 2_l a 3 



16 



32 



— a — — a 
2 8 



5 

— a 3 ) cos(kb) 



^ 2 i ' „3 



^ a' + — a 3 ) cos(2kb) + —a 3 cos(3fc6) \ . (61) 



A numerical example of this exclusive structure factor 
is presented in Fig. 8 for the case of (VO)2P20y, which 
will be discussed in Sec. VIII. 



C. Two-magnon bound state exclusive S(k) 

To evaluate the exclusive neutron scattering structure 
factor to the S = 1 two-magnon bound state we require 
the spin matrix element Eq.(|54t). This bound state is 



IV>i(*)> 



-, N d 



1(1,1) 



m, m+1/ 



+ 0(a°) |(l,l) m ,„) terms, with \m-n\ > 2, 



a (\ 
2^2^ 



-ikb 



(64) 



In the matrix element (ipi(k)\S+(x)\ipa) there is a 
— * 1 exciton coupling of the bare ground state 
to the 0(a) one-exciton component of the perturbed 
bound state, which has the same form as the leading one- 
magnon matrix element. The second contribution is the 
2^2 exciton coupling of the 0(a) |(0,0)) two-exciton 
part of the ground state to the two-exciton bare bound 
state. 

These matrix elements have quite different k- 
dependences because the — > 1 term is dimer-spin 
changing and the 2 — > 2 term is dimer-spin conserv- 
ing. From Eq.([55|) one can see that this leads to pref- 
actors of s'm(kd/2) and cos(kd/2) respectively. Since the 
one-magnon S(k) has an overall sin(fcc?/2) dependence, 
it may be possible to distinguish one-magnon and bound 
state modes by the one-magnon zeros, which occur at at 
k = 2nir/d. 

We cannot at present give a closed form result for 
the perturbative S+ bound state matrix elements for 
general k because the O(a ) mixing problem between 
states in the degenerate two-exciton manifold (discussed 
in Sec. VI) must first be solved. 

We also encountered difficulties in numerical studies 
of the bound state. We found very slow convergence of 
our projection method to this state, perhaps due to the 
presence of a nearby bulk-limit continuum, so the use 
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of our multiprecision techniques was not practical. This 
method requires implementation of an alternative itera- 
tion scheme with better convergence. Work along these 
lines is in progress; we plan to present results for bound 
state energies and matrix elements in a future publica- 
tion. 



VIII. COMPARISON WITH (VO) 2 P 2 7 

To illustrate the utility of our results we will now ap- 
ply them to the one-magnon dispersion relation ui{k) and 
exclusive neutron scattering intensity S(k) observed in 
(VO) 2 P 2 7 , "VOPO". This material is an alternating 
chain with magnetic V 4+ ion spacings of 3.2 A and 5.1 A 
along the chain pathways fl. The chains run along the 
crystallographic b direction, and there is a weaker inter- 
chain coupling J a along the "ladder" direction a, which 
gives a leading-order contribution of J a cos(fc a ) to uj{k). 

A fit of the low- lying one-magnon branch u>([0, fc& = 
k,0]) observed in VOPO to the fifth order formula 
Eq.@|9|) is shown in Fig.7. (Since VOPO has J a w 
—0.73 meV we have added this J a to the theoretical re- 
sult Eq.(31) to obtain uj[k a = 0,kh = k, k c = 0].) A 
least-squares fit gives an exchange of J = 11.0(1) meV 
and alternation of a — 0.796(4), which provides an excel- 



lent account of the data. We can similarly fit the more 
accurate ninth-order formulae for E gap and Ezb (Eqs.( 
p9| , p0| )), including the J a shift, to the measured values of 
3.1(1) meV and 15.4(3) meV; this gives consistent values 
of J = 10.92 meV and a = 0.798. 

The intensity of the modes seen in neutron scatter- 
ing is quite sensitive to the spatial geometry of the al- 
ternating chain. For VOPO we can use this relation to 
deduce which is the stronger of the two inequivalent ex- 
change paths in the chain. Fig. 8 shows the theoretical 
one-magnon intensity variation over a wide range of k. 
The dashed line is the (1 — cos(fcd)) scattering intensity 
for isolated dimers, and the solid line shows the 0(a 3 ) 
alternating-chain result Eq.(|5^) for S(k), with a = 0.8 
and d/b = 5.1 A/ 8.3 A. This d/b ratio assumes that 
the long V-PO4-V dimer, which has d — 5.1 A, has the 
stronger interaction J. The characteristic variation of 
iS(fc), which is measured in neutron scattering, allows a 
direct check of this bond assignment and is a stringent 
test of the alternating chain model as applied to VOPO 
(Ref. pi), in preparation). 



IX. SUMMARY AND CONCLUSIONS 

In summary, we have presented many new analytical 
and numerical results for the low-lying excitations of the 
spin-1/2 alternating Heisenberg chain. We introduced 
the model in Sec. I and II, and Sec. Ill reviewed previous 
work. Sec. IV. A introduced perturbation theory about 



the dimer limit, and in Sec. IV. B we used this approach 
to confirm previous results for the one-magnon dispersion 
oj{k) to third order, and extended the ground state en- 
ergy calculations to fourth order. Sec. IV. C summarized 
the expected critical behavior of the gap and ground- 
state energy as we approach the uniform chain limit. 
Relations between energies and their derivatives at the 
critical point were also derived. To complement these 
analytical calculations, in Sec.V.A we undertook numer- 
ical calculations of bulk limit energies of the ground state, 
gap, zone boundary, and 5 = and S = 1 two magnon 
bound states, using Lanczos methods on systems up to 
L = 28. As a major part of this work, in Sec.V.B we 
introduced a novel technique in which multiple precision 
calculations with a very small dimer coupling are used to 
infer perturbation expansions to high order. This con- 
firmed third and fourth order energy expansions and al- 
lowed continuation of the eo, E gap and Ezb series to 
0(a 9 ). The one-magnon dispersion relation uj(k) was 
then derived to 0(a 5 ) using this method. The expected 
critical behavior of the ground state energy and energy 
gap was compared with our numerical results in Sec.V.C. 
The ground state energy was found to be in approximate 
agreement with the expected (to within logarithmic cor- 
rections) 8 4 / 3 power law. The singlet-triplet gap however 
appears to support a power law of 5 3 ^ 4 rather than S 2 ^ 3 
over the accessible range of 5. Sec. VI considered two- 
magnon bound states. At the zone boundary k = ir/b an 
attractive potential between magnons was found to give 
rise to S = and 5 = 1 two-magnon bound states for 
all alternations. Numerical results and 0(a 2 ) analytical 
forms were given for the magnon-magnon binding energy. 
The 5 = bound state was found for all k considered, 
and numerical results for its binding energy at k = were 
also given. The 5 = 1 bound state however was found 
to exist only for a range of k near 7r/&. In Sec. VII we in- 
troduced an exclusive neutron scattering structure factor 
5(fc); this was evaluated analytically for the one-magnon 
band using the multiple precision technique and dimer 
perturbation theory The one-magnon exclusive struc- 
ture factor was determined to 0(a 3 ), and was presented 
as a modulation times the familiar isolated-dimer form 
(1 — cos(fcd)). Finally, in Sec. VIII we gave an illustrative 
application of these results to recent neutron scattering 
data on (VO)2P20y, which is dominantly an alternating 
chain. We showed that the fifth-order dispersion relation 
gives an excellent fit to the observed one-magnon band, 
and noted that the predicted exclusive structure factor 
S(k) will allow a detailed test of the alternating chain 
model as applied to (VO)2P207. 

In conclusion, we have used perturbation theory and 
numerical methods to calculate the properties of the 
ground state and one- and two-magnon states in the al- 
ternating Heisenberg chain. Using a new technique based 
on multiple precision programming we have derived high- 
order series expansions for energies and matrix elements 
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in the alternating chain. This technique is quite general 
and should be applicable to many other quasi-lD quan- 
tum spin systems. 
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TABLE I. Representative alternating chain materials. The 
parameters J = (1 + a) J/2 and 5 — (1 — a)/(l + a) are also 
commonly used (see text). 



Material J(meV) a reference 



Sri4Cll2404i 


11.2 


-0.10(1) 


25 




Cu(N0 3 ) 2 -2.5 H 2 


0.45 


0.27 






(VO) 2 P 2 7 


?a 10 


0.8 


i 




CuW0 4 


m 12 


« 0.9 






a'-NaV 2 5 


40 


0.9 


i 




CuGeOaQ 


13 


0.94 
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* There may also be important second nearest neighbor in- 
teractions in this material. 



TABLE II. Bulk limit alternating chain energies (J = 1) and two-magnon bound state binding energies, Eq.(44), extrapolated 
from L = 20, 24, 28. The change observed in going from an L = 16, 20, 24 extrapolation to L = 20, 24, 28 (an estimate of the 
systematic error) is given in parenthesis. 



a 


E /L 


U){k = 0) = Egap 


Uj(k — -K IV) 


E% l '(k = 0) 


7 — 1 — 0/7 /7\ 

E% u (k = ■n/b) 


T — i .S 1 — 1/7 / J \ 

E% (fc = 7r/&) 


0.0 [exact] 


-0.375 = -3/8 


1 


1 











0.1 


-0.375480805 


0.946279339 


1.051248884 


0.0002 


0.0456 


0.0210 


0.2 


-0.376974494 


0.885209996 


1.104980718 


0.0009(-2) 


0.0824 


0.0343 


0.3 


-0.379566321 


0.816844275(+1) 


1.161143536 


0.0025(-l) 


0.1104 


0.0403 


0.4 


t~\ o o o o f r' o f r\ 

-0.383356250 


0.74106141 (+3) 


1.219628893(+1) 


0.0049(-3) 


0.1294 


0.0402 


0.5 


-0.388465614 


0.6574777(+5) 


1.280237618(+9) 


0.0077(+34) 


0.1394 


0.0353 


0.6 


-0.395048423(-3) 


0.565296(+7) 


1.3426173(+2) 




0.1405 


0.0280(+l) 


0.7 


-0.40331243(-5) 


0.46298(+5) 


1.406138(+3) 




0.1327(-1) 


0.0207(+17) 


0.8 


-0.4135644(-8) 


0.3474(+3) 


1.46959(+4) 




0.1151(-1) 




0.9 


-0.426330(-16) 


0.2098(+17) 


1.5298(+8) 




0.0829(+2) 




1.0 


-0.44314718... 





1.5707963... 












= l/4-ln(2) 




= 7f/2 
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Figure Captions 



Figure 1. The geometry of a ID alternating chain. The internal dimer spin-spin interaction (solid line) is J\ = J 
and the dimer extent is d. The spin-spin coupling between dimers (dashed line) is J2 = aJ, and the spacing between 
dimer centers, which is the length of the unit cell, is b. A spatially uniform chain has a smaller unit cell length of 
d = b/2, which is normally called a. 

Figure 2. Ground state energy per spin eo(a) = Eq/LJ of the alternating chain. The dashed line is third-order 
perturbation theory, the solid line is ninth order, and the points are bulk limit extrapolations of Lanczos data. 

Figure 3. Singlet-triplet energy gap E gap /J of the alternating chain, as in Fig. 2. 

Figure 4. Dispersion uj{k) of the one-magnon band in the alternating chain for a = 0.2, 0.4, 0.6, 0.8 and 1.0, using 
the fifth-order dispersion relation Eqs.(j31 39). The a = 1 curve is the exact result tt\ sin(fc6/2)|/2. 



Figure 5. Critical behavior of the alternating chain as a function of 5, with J fixed. Fig. 5a shows the energy per 
spin relative to the bulk limit, eo(<5) — eo(0), and a fit to ciS 4 ^ 3 which gives c\ = 0.3134. Fig. 5b shows E gav j J and a 
fit to ci£ 3 / 4 which gives ci = 2.003. 

Figure 6. Binding energies of the S — and S — 1 two-magnon bound states at k = n/b (band maximum). The 
points are bulk-limit extrapolations of Lanczos data (Table II) and the lines are second order perturbation theory, 
Eq.©. 



Figure 7. A fit of the alternating chain dispersion relation Eqs.(pT 39) to the one-magnon dispersion observed 
in neutron scattering from (VO)2P207 [Q, as discussed in Sec. VII. The fitted parameters are J = 11.0 meV and 
a = 0.796. 

Figure 8. The predicted neutron scattering intensity (exclusive structure factor) S(k) for the one-magnon mode in 
(VO)2P207, Eq.(|6l|). The parameters are a = 0.8 (from a fit to the dispersion), d = 5.1 A and b = 8.3 A. The dashed 
line shows the isolated dimer form (1 — cos(fcfi)) for comparison. 
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1 a) real space lattice 
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